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Abstract 

Proton acceleration by ultra-intense laser pulse irradiating a target with cross-section smaller 
than the laser spot size and connected to a parabolic density channel is investigated. The target 
splits the laser into two parallel propagating parts, which snowplow the back-side plasma electrons 
along their paths, creating two adjacent parallel wakes and an intense return current in the gap 
between them. The radiation-pressure pre-accelerated target protons trapped in the wake fields now 
undergo acceleration as well as collimation by the quasistatic wake electrostatic and magnetic fields. 
Particle-in-cell (PIC) simulation shows that stable long-distance acceleration can be realized, and 
a 30 fs monoenergetic ion beam of > 10 GeV peak energy and < 2.0° divergence can be produced 
by a 9.8 x 10^-*^ W/cm^ circularly polarized laser pulse. 

PACS numbers: 52.38.Kd, 41.75.Jv, 52.35.Mw,52.59.-f 
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Laser-driven ion acceleration has been of much recent interest because of its broad scien- 
tific and technical applications. Protons with tens to hundreds MeV energies are useful for 
cancer therapy [1], high resolution imaging , fast ignition in inertial confinement fusion 
PI, etc., and ions with still higher energies are relevant to high energy physics research l^l, 
laboratory modeling of astrophysical phenomena jo^, etc. 

Radiation pressure acceleration (RPA) is promising for obtaining high-quality ion beams 
efl&ciently I?- 9]. Existing studies have shown that GeV proton beams can be obtained 
by RPA using CP lasers with intensity above 10^^ W/cm^ I id. However, because of 
transverse instabilities and hole-boring by the laser pulse (l0l-ll2l| , the acceleration length is 
rather limited and it is difficult to enhance proton energy without still higher laser inten- 
sity. Recently, it has been shown that tens GeV proton beams Q be generated 



by a moving double-layer in a possibly unlimited acceleration regime, or in a two-phase 
acceleration regime Q by ultra-relativistic lasers. In classical RPA and laser wakefield 
acceleration (LWFA), however, the accelerated ion bunch tends to diverge because of the 
ubiquitous space-charge field (l7l-ll9[|. The resulting defocusing leads to very low ion fiux. 
Moreover, the effective laser-plasma interaction distance is also limited by laser diffraction. 
As a result, ions cannot remain collimated and be stably accelerated over a long distance. 

In this Letter, we propose a stable proton acceleration scheme using an ultra-intense 
circularly polarized (CP) laser pulse and a mass-limited target (MLT) [20] connected to a 
parabolic underdense plasma channel. The target splits the laser pulse into two separate 
but adjacent parts, which are then refract ively guided by the preformed underdense plasma 
channel behind it. The quasistatic longitudinal accelerating and transverse coUimating fields 
associated with the electron return current at the center of the channel can stably trap, 
coUimate, and further accelerate the RPA target protons by LWFA over a long distance. In 
contrast to the existing acceleration schemes, both proton defocusing and laser diffraction 
do not occur. PIC simulations show that ultrashort (30 fs) monoenergetic proton beams 
with peak energy > 10 GeV and divergence angle < 2.0° can be obtained with a ^ 10^^ 
W/cm^ circularly polarized (CP) laser pulse. 

To demonstrate the proposed scheme we use the two-dimensional (2D) PIC code KLAP 
fiol . I21I . The CP laser pulse is of wavelength A = 1 /xm and peak intensity Iq = 9.8 x 10^^ 
W/cm^ (or the normalized laser parameter ag = 60). It has a super-Gaussian radial profile 
with spot size 2a = 40A and a 18T fiattop envelope, where T is the laser period. The fully 
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FIG. 1: (Color online) (a) Laser electric field ^ E'^ + E^^ normalized by £{). (b) Electron density 
normalized by Uc. (c) Current density jx normalized by ericC. (d) Spatial distribution of transverse 
quasistatic field Eys — B^s (in the ±y directions), at t ~ 500T. 

ionized uniform hydrogen target of density of = 20nc, where ric is the critical density, is 
located in 50A < x < 50. 8A, — lOA < y < lOA. Behind the hydrogen target is an underdense 
plasma channel with a parabolic density profile n = rii + /\n{y'^ /y^)^ where rii = 0.15nc, 
An = 2nc, yo = 50A, and it is located in 50. 8A < x < 600A, — 40A < y < 40A. The simulation 
box is of size lOOA x lOOA, corresponding to a moving window consisting of 3200 x 1000 
cells, with each cell containing 800 macroparticles for the target and 6 for the underdense 
channel. 

Figs. [11(a) and (b) show the laser intensity and electron density, respectively, at t ^ 500T. 
The transverse size h of the target is smaller than the laser spot size a. The incoming laser 
pulse is split into two by the target. After passing the latter they continue to propagate 
forward in the preformed backside underdense plasma channel as adjacent twin pulses. The 
ponderomotive force of the latter expels the local plasma electrons transversely outward. 
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creating two adjacent wakes, or twin wakes, which trap the RPA target protons. Fig. [T](c) 
shows the longitudinal current density jx at t ^ 500T. One can see that there is an intense 
electron return-current sheet between the two wakes. Since the longitudinal velocity Vx of the 
trapped protons is near the light speed c, the transverse field experienced by the energetic 
protons can be written as Scy = EygCy + c~^VxexC x B^sCz ^ {Eyg — Bzs)^y^ where Eys and 
Bzs are the self-generated electrostatic and magnetic fields, respectively. Fig. [l](d) shows 
the total transverse quasistatic field £ = Eys — Bzs normalized by = vn^ujcje^ where mg 
and — e are the electron mass and charge, and uj and c are the laser frequency and light 
speed, at t ^ 500T. One can see that the twin wakes act as a moving potential well for 
confining, stabilizing, and further accelerating (now by LWFA) the RPA protons for a long 
distance in underdense plasma. In addition, since the target electrons also move with the 
twin wakes, the expelled plasma electrons are prevented from rapidly refilling the back side 
of the channel, so that the latter remains open for a long distance. This scenario diflFers 
significantly from that of the classical bubble/ blowout acceleration regimes [22]. 

When the initial laser pulse is split by the target, the eflFective spot size g\ of each of the 
resulting twin light pulses is less than that of the original pulse. Similarly, the divergence 
angle = (Ji/Zr is larger and the Rayleigh length Zr = Traf/X smaller. Earlier results 
have also shown that relativistic self-focusing cannot prevent diffraction of intense light 
pulses 3, Q. Thus, without the underdense plasma channel at the back of the target the 
laser-plasma interaction distance would be very limited. 

However, this undesirable scenario is completely changed by the presence of the preformed 
parabolic density channel, here with the profile n = [rii + /\n{y^ /y^)]^ at the back of the 
target. The dispersion relation for relativistic light waves is uj^ = c^k'^ + cl;^/7, where k 
and (joI/j are the wave number and effective background plasma frequency, respectively, and 



7 ^ \/T+~a? is the relativistic factor associated with the electron quiver motion. For an 
ultra-intense CP laser pulse (a 1), the index of refraction is 
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so that refractive light guiding (RLG) occurs if dyfj < 0. Light diffraction can then be 
canceled by RLG if the maximum beam focusing angle is larger than the divergence angle 
Od, or 
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FIG. 2: (Color online) Density of the accelerated protons in units of ric at (a) t ^ lOOT and (c) 
t ^ 500r. (b) Angular distribution of the protons at t ^ 600T. The dashed lines correspond to 
the divergence angle ±2.5°. (d) Evolution of the proton energy spectrum. 

where the three terms of the inequality represent relativistic self-focusing guiding, preformed 
density channel RLG, and light diffraction, respectively. We note that this relation differs 
considerably from that of the a <C 1 case. For the parameters under consideration, the self- 
focusing and RLG terms in Eq. ([2]) have the values ^ 0.0021 and 0.013, respectively. The 
latter alone is already much larger than the diffraction term, which has the value 0.0041. One 
can also see that without RLG the proton bunch would diverge. Fig. [l](a) shows the profile 
of the light field at t ^ 500T. We see that RLG by the preformed parabolic density channel 
indeed results in long distance propagation of the twin light pulses and the corresponding 
twin wakes. In fact, the laser-plasma interaction distance is effectively increased to ^ 7Zr. 
Moreover, the focusing effect of RLG also doubles the intensity of the twin pulses to ai = 120, 
as can be seen in Fig. [I](a). 

We now consider the evolution of the pre-accelerated protons in more detail. From Figs. 
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FIG. 3: (Color online) Evolution of the (a) maximum proton energy, and (b) number of protons 
in the axis region(|y| < 3A). 

[2](a) and (c) one sees that the relativistic protons are confined along the original laser axis and 
they can propagate a long distance without divergence. Fig. [T](d) shows that the maximum 
quasistatic transverse (or focusing) field (?max induced by the charge separation and the 
return current at the center of the channel is about GSq^ agreeing well with the theoretical 



prediction S^s,^/Sq ^ (7r/2)(aini/nc)'^/^[22|. The angular distribution of the accelerated 
protons at t ^ 600T is given in Fig. [2](b), which shows that the average emission angle for 
protons with energy greater than 8 GeV is less than 2.0°. That is, the output proton beam 
is well coUimated by the transverse quasistatic field. 

Fig. [3](a) shows the evolution of the maximum proton energy. One can see the two 
stages of the acceleration process, namely RPA and LWFA. In the initial RPA regime, the 
maximum proton velocity reaches 0.945c, or 1.93 GeV, at t ^ lOOT. The pre-accelerated 
protons then undergo LWFA in the twin wakes until the light pulses are almost depleted at 
t ^ 600T. At that time, the trapped proton bunch approaches the front of the twin wakes, 
indicating that the dephasing length is roughly the same as the pump depletion length. We 
also see that the maximum proton energy increases rapidly with time and then saturates at 
t ^ 600T. Fig. [2](d) shows the evolution of energy spectrum. One can see that the latter 
improves with time and eventually achieves a sharp peak at about 12 GeV, which is also 



consistent with that of a theoretical model 



15| Wmax - {alnc/ni)meC^/6. Fig. [3](b) shows 



that the number of protons in the region \y\ < 3X first decreases and then becomes constant 
at 1.5 X 10^^, indicating that the transverse spreading effect for the energetic protons is well 
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FIG. 4: (Color online) Electron charge distributions at t ^ GOOT for the (a) single-wake (when the 
condition /i < a is not satisfied) and (c) twin-wakes cases. The corresponding angular distributions 
of the protons are shown in (b) and (d), respectively. 

controlled in the backside channel. Their total energy is above 500 J and the laser conversion 
efficiency is about 10%. The accelerated protons are well compressed in the phase space and 
the quasi-monoenergetic pulse has a duration of about 30 fs. Such a short proton pulse can 
be used to excite a plasma wakefield that can in turn accelerate the plasma electrons [23] . 

A key issue in the proposed scheme is the formation of a stable proton focusing and 
collimating twin- wakes structure that moves at relativistic speed with the target-modulated 
laser light. When the transverse size of target is too large comparing to the laser spot size, 
the initial laser pulse can fully interact with the target. The laser is not split into two 
parts and no twin wakes are formed. Without the transverse confinement in the underdense 
region, the RPA protons would diverge. This can be seen in the simulation by setting 
h = 2a = 40/xm, with the other laser and plasma parameters the same as in the preceding 
simulations. Fig. [4](a) shows that instead of the stable twin wakes (see Fig. [11(b)) a single 
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channel is generated by the laser pulse. Due to the defocusing effect of the space-charge field, 
the RPA protons quickly diverge in the wake, as can be seen in Fig. [11(b) for the emission 
angle (almost 20°) at t ^ 600T. In addition, protons that deviate from the axis region do 
not experience LWFA. The maximum energy of the accelerated protons is about 5 GeV, 
which is lower than that in the the stable twin- wakes case (see Fig. [21(d)). The quality of the 
final proton bunch is thus relatively poor in terms of energy gain as well as emission angle. 

The condition ([2j) for guiding the laser light in the underdense plasma is also important. If 
the maximum focusing angle is less than the beam divergence angle, RLG of the twin pulses 
in the density channel cannot occur. Instead, the acceleration process becomes unstable 

nn 

(13l . Il4j |. To see this, we set An = (uniform plasma behind the target), with the other 
parameters unchanged. That is, the density- channel guiding effect is removed. Since the 
relativistic self-focusing term is less than the beam diffraction term. Fig. [11(c) shows the 
density distribution of electrons in this case. We see that although a twin- wakes structure is 
still formed at the beginning, it does not survive for a long time. Only a fraction of the RPA 
protons is further accelerated and for a much shorter distance than that in the preformed 
channel. In addition, the average divergence angle is now rather large, namely 14°, as shown 
in Fig.[lKd). 

In conclusion, an efficient laser-driven proton acceleration scheme using a small target 
with a preformed parabolic underdense channel behind is proposed. The laser pulse drives 
RPA protons and is split by the small target into two parts, which excite behind them a 
twin-wakes structure in the channel. The pre-accelerated RPA target protons are efficiently 
collimated and further accelerated by LWFA in the wake. With this scheme, both light 
diffraction and proton defocusing are overcome, making long-distance acceleration possible. 
The conditions and characteristics of proton acceleration from the PIC simulations agree well 
with that from the analytical estimates. It is shown that ultrashort (30 fs) highly-collimated 
monoenergetic proton bunches with peak energy > 10 GeV can be produced with a 60 fs CP 
laser pulse of intensity 9.8 x lO^-*^ W/cm^. These bunches may be superior in fs metrology 
in determining the validity of the OPERA type experiments for neutrino velocity (2^. Such 
a short and collimated proton pulse can also be promising to excite a plasma wakefield for 
TeV electron acceleration. 
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